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Abstract 



The result of a search at the LHC for heavy stable charged particles produced in pp 
collisions at \/s = 7 TeV is described. The data sample was collected with the CMS 
detector and corresponds to an integrated luminosity of 3.1 pb _1 . Momentum and 
ionization-energy-loss measurements in the inner tracker detector are used to iden- 
tify tracks compatible with heavy slow-moving particles. Additionally tracks passing 
muon identification requirements are also analyzed for the same signature. In each 
case, no candidate passes the selection, with an expected background of less than 
0.1 events. A lower limit at the 95% confidence level on the mass of a stable gluino 
is set at 398 GeV/ c 2 , using a conventional model of nuclear interactions that allows 
charged hadrons containing this particle to reach the muon detectors. A lower limit 
of 311 GeV/c 2 is also set for a stable gluino in a conservative scenario of complete 
charge suppression, where any hadron containing this particle becomes neutral be- 
fore reaching the muon detectors. 



Submitted to the journal of High Energy Physics 



*See Appendix A for the list of collaboration members 



1 



1 Introduction 

Heavy stable (or quasi-stable) charged particles (HSCPs) appear in various extensions of the 
standard model (SM) HHSl. If the lifetime of an HSCP produced at the Large Hadron Collider 
(LHC) is longer than a few nanoseconds, the particle will travel over distances that are compa- 
rable or larger than the size of a typical particle detector. In addition, if the HSCP mass is >100 
GeV/ c 2 , a significant fraction of these particles will have a velocity, j6 = v/c, smaller than 0.9. 
These HSCPs will be directly observable through the distinctive signature of a high momentum 
(p) particle with an anomalously large rate of energy loss through ionization (dE/ dx). 

Previous collider searches for HSCPs have often been performed under the assumption that 
these particles lose energy primarily through low-momentum-transfer interactions, even if they 
are strongly interacting, and are therefore likely to reach the outer muon systems of the detec- 
tors and be identified as muons ||9TH2|. The interactions with matter experienced by a strongly- 
interacting HSCP, which is expected to form a bound state (R-hadron) [13| in the process of 
hadronization, can lead to it flipping the sign of its electric charge or becoming neutral. A 
recent study [14J on the modeling of nuclear interactions of HSCPs traveling through matter, 
favours a scenario of charge suppression. In this model the probability is close to unity for 
an R-hadron containing a gluino, g (the supersymmetric partner of the gluon), or a supersym- 
metric bottom squark, to emerge as a neutral particle after traversing an amount of material 
typical of the detectors operating at LEP, the Tevatron, or LHC. If this model is correct, the ma- 
jority of these HSCPs would not be observed in the muon system of a typical collider detector. 
Experimental strategies that do not rely on the muon-like behavior for the HSCPs are there- 
fore of great importance. For instance, searches have been performed for very slow (ft < 0.4) 
R-hadrons containing a gluino brought to rest in the detector llT5lll"6l . 

In this article we present a search for HSCPs produced in pp collisions at yjs = 7 TeV at the 
LHC with the Compact Muon Solenoid (CMS) detector [17]. The search is based on the data 
sample collected between April and August 2010 corresponding to an integrated luminosity 
of 3.06 pb _1 . We use triggers requiring: a high-transverse-momentum muon (pj > 9 GeV/c); 
or a dimuon pair (pj > 3 GeV/c for each muon); or calorimeter-based missing transverse 
energy (E™ 1SS > 100 GeV), to search for HSCPs failing muon identification or emerging mainly 
as neutral particles after traversing the calorimeters; or a high-transverse-energy jet (Ej > 100 
GeV) to search for HSCPs accompanied by substantial hadronic activity. The analysis makes 
use of two approaches. In a first selection, referred to as "tracker-only" , the HSCP candidates are 
searched for as individual tracks reconstructed in the inner tracker detector with large dE/ dx 
and pj. A second selection, referred to as "tracker-plus-muon" , additionally requires that the 
track is identified as a muon in the outer muon detector. For both selections, the mass of the 
candidate is calculated from the measured p and dE/ dx. 

2 The CMS Detector 

The central feature of the CMS detector is a 3.8 T superconducting solenoid of 6 m internal di- 
ameter surrounding a silicon pixel and strip tracker, a crystal electromagnetic calorimeter, and 
a brass-scintillator hadronic calorimeter. Muons are measured in gaseous detectors embedded 
in the iron return yoke. Centrally produced charged particles are measured in the tracker by 
three layers of silicon pixel detectors, followed by ten microstrip layers. At pseudorapidities 
(t] = — lntan(0/2), where 9 is the polar angle measured with respect to the beam direction) 
above ~ 1.5, particles are tracked in two pixel and twelve strip layers arranged in disks per- 
pendicular to the beam axis. In this analysis, the dE/dx measurement is based only on the 
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3 Candidate Selection and Background Estimation 



information from the silicon strip detectors. The dE/ dx measurement precision is limited by 
the silicon strip analogue-to-digital converter (ADC) modules that are characterized by a max- 
imum number of counts per channel corresponding to about three times the average charge 
released by a minimum-ionizing particle (MIP) in 300 }im of silicon. This is the thickness of 
the modules mounted in the innermost silicon strip central layers. The pj resolution for tracks 
measured in the central (forward) region of the silicon tracker is 1% (2%) for pj values up to 
50 GeV/ c and degrades to 10% (20%) at pj values of 1 TeV/ c. The trigger and reconstruction 
efficiencies for HSCPs in the muon detectors are limited by the requirements on the arrival time 
of the particles at the muon system. These requirements affect the efficiency for detecting slow 
HSCPs. The dependence of the muon trigger efficiency on the particle velocity (ft) is studied 
using data and Monte Carlo (MC) simulations and found to decrease, below f> = 0.7. The 
muon trigger becomes completely inefficient at f> = 0.5. A much more detailed description of 
the CMS apparatus can be found elsewhere |T7|. 



Candidate HSCPs are pre-selected by requiring a track with \f]\ < 2.5, pj > 15 GeV/ c, relative 
uncertainty on the pj less than 15%, and transverse (longitudinal) impact parameter with re- 
spect to the reconstructed primary collision vertex less than 0.25 (2.0) cm. Candidate tracks are 
also required to have at least three measurements in the silicon-strip detector. For the tracker- 
plus-muon selection, we additionally require the track to be compatible with track segments re- 
constructed in the muon system. As an estimator of the degree of compatibility of the observed 
charge measurements with the MIP hypothesis, a modified version of the Smirnov-Cramer- 
von Mises (181 H9| discriminant is used (the modification applied to the original form of the 
discriminant eliminates the sensitivity to incompatibility with the MIP hypothesis due to low 
ionization): 



where N is the number of charge measurements in the silicon-strip detectors, P, is the probabil- 
ity for a MIP to produce a charge smaller or equal to the z-th charge measurement for the ob- 
served path length in the detector, and the sum is over the track measurements ordered in terms 
of increasing P,. The charge probability density function used to calculate P, is obtained using 
tracks with p > 5 GeV/ c in events collected with a minimum bias trigger. Non-relativistic 
HSCP candidates will have the value of the discriminant I as approaching unity. Figure [l] shows 
normalized distributions of pj and I as in data and two MC samples, for candidates passing 
the tracker-only pre-selection. The first MC sample contains events from QCD processes. The 
second MC sample contains signal events from pair-production of stable g with a mass of 200 
GeV/ c 2 . Both samples are generated with the PYTHIA v6.422 [20] MC package. More details on 
the simulation of the signal sample will be given below. The MC QCD simulations are found 
to reproduce the data, and the simulated signal is clearly separated. Because of the limited 
number of available simulated events with low transverse-momentum transfers, the MC QCD 
distributions display bin-to-bin variations in the size of the statistical errors. 

The most probable value of the particle dE/ dx is determined using a harmonic estimator I/, of 
grade k = —2: 



where c, is the charge per unit path length in the detector of the z'-th measurement for a given 
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Figure 1: Normalized distributions of pj (left) and 7 flS (right) in data and two MC samples, for 
candidates passing the tracker-only pre-selection. The two MC samples contain events from 
QCD processes and from pair-production of g with a mass of 200 GeV/ c 2 , respectively. 



track. In order to estimate the mass (m) of highly ionizing particles, the following relationship 
between If,, p, and m is assumed: 

2 

m 

I h = K^ + C. (3) 

Equation [3] reproduces the Bethe-Bloch formula |2l| with an accuracy of better than 1% in the 
range 0.4 < /3 < 0.9, which corresponds to 1.1 < (dE/dx)/ (dE / dx) mip < 4.0. The empirical 
parameters K and C are determined from data using a sample of low-momentum protons, for 
which the fitted values are K = 2.579 ± 0.001 MeV cm" 1 c 2 and C = 2.557 ± 0.001 MeV cm" 1 , 
and the mass resolution is 7%. The reconstructed mass distribution for kaons and protons is in 
very good agreement with the one obtained from MC following this procedure |22] . For masses 
above 100 GeV/c 2 , the mass resolution is expected to worsen because of the deterioration of the 
momentum resolution and because of the limit on the maximum charge that can be measured 
by the silicon strip tracker ADCs, which also affects the mass scale. For a 300 GeV/c 2 HSCP, 
the mass resolution is 12% and the reconstructed peak position is at 265 GeV/ c 2 . 

The search is performed as a counting experiment. Signal candidates are required to have I as 
and Pt greater than threshold values and the mass to be in the range of 75 to 2000 GeV/c 2 , 
allowing sensitivity to HSCP masses as low as 100 GeV/c 2 . The I as distribution for the pre- 
selected tracks, and in particular its tail, depends strongly on the number of charge measure- 
ments on the track. Thus, to increase the sensitivity of the search, pre-selected tracks are di- 
vided into subsamples according to the number of silicon strip measurements. The I as (pj) 
threshold in each subsample is determined by requiring a constant efficiency on data for all 
subsamples, when the threshold is applied separately. A method that exploits the absence 
of correlation between the pj and dE/ dx measurements in data is used to estimate the back- 
ground from MIPs. In a given subsample j, the number of tracks that are expected to pass both 
the final pj and l as thresholds set for the subsample is estimated as D, = B,C,/ At, where Aj is 
the number of tracks that fail both the I as and pj selections and Bj (Cj) is the number of tracks 
that pass only the I as (pj) selection. The By and Cj tracks are then used to form a binned proba- 
bility density function in If, (p) for the Dy tracks. Finally, using the mass determination (Eq.|3jl, 
the full mass spectrum of the background in the signal region D is predicted. 
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4 Results 



By comparing the predicted and observed number of tracks for several very loose selections in 
a control region of the mass spectrum, corresponding to masses below 75 GeV/ c 2 , the predic- 
tion is found to underestimate systematically the observation by 12% (5%) for the tracker-only 
(tracker-plus-muon) selection. After correcting the predicted background by this amount, the 
remaining background systematic uncertainty is conservatively estimated as twice the r.m.s. 
of the prediction-to-observation ratio distribution The resulting uncertainty on the predicted 
background is 14% (17%). 

As significant background rejection can be obtained without a sizable effect on the signal effi- 
ciency, the final selection is optimized by requiring the total expected background in the search 
region to be ~ 0.05 events. This low-background choice optimizes the discovery potential even 
if just a handful of events are observed, and at the same time maintains significant exclusion 
sensitivity in the case that no events are observed. 



4 Results 

In addition to the final "tight" selection, the result of a "loose" selection is reported in Table [l] 
The loose selection retains a relatively large number of background candidates and allows us 
to compare the background prediction with the observed data. Figure|2]shows good agreement 
between the observed and predicted mass spectrum obtained using the loose selection for the 
tracker-plus-muon and tracker-only candidates. 




Mass (GeV/c 2 ) Mass (GeV/c 2 ) 



Figure 2: Mass spectrum for the loose selection defined in Table |l| for the tracker-plus-muon 
(left) and tracker-only (right) candidates. Shown are: observed spectrum (black dots with 
the error bars), data-based predicted background spectrum (red triangles) with its uncertainty 
(green band) and the spectrum predicted by MC for a signal of pair-produced stable g with a 
mass of 400 (left) and 300 (right) GeV/c 2 (blue histogram). 



The results of the search with the final selection are also presented in Table [T] No candidate 
HSCP track is observed in either the tracker-only or tracker-plus-muon analysis. 

Given the null result, cross section upper limits at the 95% C.L. are set on the HSCP production 
for two benchmark scenarios: direct production of g pairs and ti pairs. For a given mass, the 
cross section for g production is expected to be much larger than that for ti production at both 
the Tevatron and the LHC. Thus higher mass limits can be set for the former at both machines. 
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Table 1: Selections used in the analysis and results of the search. The tracker-plus-muon and 
tracker-only selections are labeled as "Mu" and "Tk", respectively. As explained in the text, 
the actual I as (Pt) thresholds are determined in the various subsamples by the requirement of a 
constant efficiency for candidate selection, ej (e pT ). These thresholds, indicated by I'" s ' n (pj m ), 
are therefore reported as a range of values. Expected and observed number of candidates in 
the signal region are reported in the "Expected" and "Observed" rows, respectively. Top: loose 
selection. Bottom: tight selection. 



LOOSE 


Mu 


Tk 




3.2 x 1(T 2 


1.0 x 10~ 2 




0.049 - 0.162 


0.007 - 0.278 


€ Pt 


1.0 x 10 _1 


3.2 x 10~ 2 


p™ n (GeV/c) 


34-36 


59-62 


Expected 


281±2(sfat) ±49(sysf.) 


426±l(stat.)±62(syst.) 


Observed 


307 


452 


TIGHT 


Mu 


Tk 


ei 


1.0 x 10~ 4 


1.0 x 10~ 4 


rmin 
1 as 


0.184 - 0.782 


0.186 - 0.784 


€ Pt 


1.0 x 10~ 3 


3.2 x 10~ 4 


p™ n (GeV/c) 


115-118 


154 - 210 


Expected 


0.025 ± 0.002(sfflf.) ± 0.004(syst.) 


0.074 ± 0.002(stat.) ± 0.011 (syst.) 


Observed 









However, as the mass of a produced particle increases, the ratio of the production cross section 
at the LHC to that at the Tevatron increases. For g masses in the region of 350 GeV/ c 2 , the 
increase in relative cross section outweighs the difference in integrated luminosity between the 
current Tevatron and LHC data sets, enabling the LHC to set the most sensitive limits on the 
search for g. 

Events with pair production of g and t\, with mass values in the range 130-900 GeV/ c 2 , are gen- 
erated with pythia in order to compute the efficiency of our selection on these signals. The ti 
and g are treated as stable in all these samples and their hadronization is performed by pythia. 
A parameter relevant to the g pair production, and not to the ti pair production, is the fraction, 
/, of produced g hadronizing into a g-gluon state (R-gluonball). This fraction is an unknown 
parameter of the hadronization model and affects the fraction of R-hadrons that are neutral at 
production, which in turn affects the detection efficiency. In this study, results are obtained for 
two different values of/, 0.1 and 0.5, to show the effect of the hadronization model uncertainty 
on the sensitivity of the search. The interactions of the HSCPs with the CMS apparatus and the 
detector response are simulated in detail with the GEANT4 v9.2 Il23ll24| toolkit. The R-hadron 
strong interactions with matter are modeled as in Ref. fl25ll26| . This model, like a number of 
others (fl4ll27ll29l , assumes that the probability of an interaction between the heavy parton and 
a quark in the target nucleon is low since the cross section varies with the inverse square of the 
parton mass according to perturbative QCD. The adopted model chooses a pragmatic approach 
based on analogy with observed low energy hadron scattering. However, given the very large 
uncertainties on the dynamics underlying R-hadron interactions, an extremely pessimistic sce- 
nario of complete charge suppression, where each nuclear interaction suffered by the R-hadron 
causes it to become neutral, is also considered. The tracker-only selection is expected to have 
sensitivity even in such a scenario. The total signal efficiency is reported in Table [2] for some 
combinations of models and selections. Relatively small differences are found between the 
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4 Results 



Table 2: Total signal selection efficiency and cross section upper limits for different combina- 
tions of models and selections: pair production of supersymmetric stop and gluinos; tracker- 
plus-muon (Mu) and tracker-only (Tk) selections; different fractions, /, of K-gluonball states 
produced after hadronization and charge suppression (ch. suppr.) scenario. 



gluino mass (GeV/ c 2 ) 


200 


300 


400 


500 


600 


900 


Theoretical cross section (pb) 


606 


57.2 


8.98 


1.87 


0.470 


0.0130 


Mu; f=0.1 














Total efficiency (%) 


7.17 


10.4 


13.1 


15.1 


14.5 


9.18 


Expected 95% C.L. limit (pb) 


15.1 


10.4 


8.25 


7.16 


7.47 


11.8 


Observed 95% C.L. limit (pb) 


14.5 


9.98 


7.92 


6.88 


7.17 


11.3 


Mu; f=0.5; 














Total efficiency (%) 


3.84 


5.46 


7.03 


8.23 


8.10 


4.98 


Expected 95% C.L. limit (pb) 


28.2 


19.8 


15.4 


13.1 


13.3 


21.7 


Observed 95% C.L. limit (pb) 


27.1 


19.0 


14.8 


12.6 


12.8 


20.9 


Tk; f=0.1; ch. suppr. 














Total efficiency (%) 


0.59 


2.44 


4.16 


6.39 


8.60 


7.66 


Expected 95% C.L. limit (pb) 


188 


45.5 


26.7 


17.4 


12.9 


14.5 


Observed 95% C.L. limit (pb) 


176 


42.6 


25.0 


16.2 


12.1 


13.6 


stop mass (GeV/ c 2 ) 


130 


200 


300 


500 


800 




Theoretical cross section (pb) 


120 


13.0 


1.31 


0.0480 


0.00110 




Mu; 














Total efficiency (%) 


2.99 


9.50 


14.7 


19.6 


14.0 




Expected 95% C.L. limit (pb) 


36.1 


11.4 


7.35 


5.52 


7.71 




Observed 95% C.L. limit (pb) 


34.7 


10.9 


7.06 


5.30 


7.39 




Tk; ch. suppr. 














Total efficiency (%) 


0.02 


1.19 


3.55 


7.27 


7.68 




Expected 95% C.L. limit (pb) 


5540 


93.2 


31.3 


15.3 


14.5 




Observed 95% C.L. limit (pb) 


5180 


87.2 


29.2 


14.3 


13.5 
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Table 3: Sources of systematic errors and corresponding relative uncertainties. 



Source of Systematic Error 


Relative Uncertainty (%) 


Theoretical cross section 


10-25 


Integrated luminosity 


11 


Trigger efficiency 


12 


Muon reconstruction efficiency 


5 


Track reconstruction efficiency 


< 5 


Momentum scale 


< 5 


Ionization energy loss scale 


< 3 


Total uncertainty on signal acceptance 


15 



tracker-plus-muon and tracker-only selection except in the charge suppression scenario, where 
the tracker-plus-muon selection is completely inefficient. 

This analysis is found to be complementary to the search for long-lived stopped particles pre- 
sented in [16] . Indeed, for the case of g with / = 0.1 and mass values below 500 GeV/ c 2 , the 
fraction of HSCPs that have /3 < 0.4 and pass the final selection is less than 0.5%. Therefore the 
two analyses explore different ranges of produced particle velocities with no overlap. 

The main sources of systematic uncertainty affecting the results presented in the following 
are summarized in Table [3] The uncertainty on the signal selection efficiency is estimated to 
be 15% for all considered combinations of models and selections. The main source of this 
uncertainty is an assumed 10% uncertainty on the jet energy scale 11301 , which affects both the 
jet and Ej lss trigger efficiency by about 10%. In a more recent study |3Tf, the estimate of the 
uncertainty on the jet energy scale has been reduced by a factor of two. However, in this 
analysis we have conservatively chosen to retain the earlier estimate of 10%. The uncertainty 
on the muon trigger efficiency and the impact of an imperfect simulation of the synchronization 
of the muon trigger and readout electronics are studied with data and MC. They result in an 
overall uncertainty on the signal selection efficiency of less than 5%. The uncertainty on the 
offline track reconstruction efficiency |I32U33II , track momentum scale 113411 and ionization energy 
loss scale [22J is also found to yield no more than 5% uncertainty on the overall signal selection 
efficiency. The uncertainty on the absolute value of the integrated luminosity is estimated to be 
11% 135].' 

The upper limit on the cross section is computed at 95% C.L. using a Bayesian method with a 
flat signal prior and a log-normal prior used for integration over the nuisance parameters ffl8l 
[l9j|2T]|. In order to obtain a conservative upper limit we set the expected background to zero. 
The tracker-plus-muon selection provides better limits than the tracker-only for all scenarios 
but the one with complete charge suppression. For each considered scenario, the cross section 
upper limit obtained with the most sensitive selection is reported in Tableland Fig. [3j along 
with the theoretical predictions for g and ti pair production computed at next-to-leading order 
(NLO) + next- to-leading log (NLL) |36ti39l using the PROSPINO v2 program |40j. The g theo- 
retical predictions refer to the case where the squarks and gluino are degenerate in mass. In 
the heavy squark limit these cross sections are about 10% higher. For the case of fi, beyond 
LO, the cross section does not only depend on the ti mass, but also, though to a much lesser 
extent ]4T|, on the g mass, the average mass of the first and second generation squarks and the 
stop mixing angle. For this reason, the ti theoretical predictions reported in Tableland Fig. [3] 
refer to the SPSla' benchmark scenario ]42|- All systematic uncertainties discussed above are 
included in the cross section upper limits reported in Tableland Fig. |3] From the intersection 
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of the cross section limit curve and the lower edge of the theoretical cross section band we set 
a 95% C.L. lower limit of 398 (357) GeV/c 2 on the mass of pair-produced g with / = 0.1(0.5), 
using the tracker-plus-muon selection. The analogous limit on the ti mass is 202 GeV/c 2 . In 
the charge suppression scenario we set, with the tracker-only selection, a 95% C.L. g mass limit 
of311GeV/c 2 for/ = 0.1. 

5 Conclusions 

In summary, the CMS detector has been used to identify highly ionizing, high-p?- particles 
and measure their masses. Two searches have been conducted: a very inclusive and model 
independent one that uses highly-ionizing tracks reconstructed in the inner tracker detector, 
and another requiring also that these tracks be identified in the CMS muon system. In each case, 
the observed distribution of the candidate masses is consistent with the expected background. 
We have set lower limits on masses of stable strongly interacting supersymmetric particles. For 
the case of g with / = 0.1 and ti, a lower mass limit of 398 and 202 GeV/c 2 , respectively, is 
set at the 95% C.L. with the analysis that uses muon identification. In a pessimistic scenario of 
complete charge suppression the above g mass limit is reduced to 311 GeV/ c 2 and is obtained 
with the tracker-only selection. The limits presented here on stable g are the most restrictive to 
date. 
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